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ABSTRACT 

In the framework of the fluctuation-dissipation approach to dynamical friction, we 
derive an expression giving the orbital energy exchange experienced by a compound 
body as it moves interacting with a non homogeneous discrete background. The body 
is assumed to be composed of particles endowed with a velocity spectrum and with 
a non homogeneous spatial distribution. The Chandrasekhar formula is recovered in 
the limit of a point-like satellite with zero velocity dispersion and infinite temperature 
moving through an homogeneous infinite medium. In this same limit, but dropping 
the zero satellite velocity dispersion (as) condition, the orbital energy loss is found to 
be smaller than in the as — case by a factor of up to an order of magnitude in some 
situations. 
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1 INTRODUCTION 

A satellite moving through a field of gravitating particles 
experiences a dissipative frictional force known as dynam- 
ical friction. It can be understood in terms of the satellite 
wake, that exerts a drag force on the satellite itself (Kalnajs 
1972; Binney & Tremaine 1987; Weinberg 1989; Bekenstein 
& Zamir 1991). It can also been understood in terms of 
the underlying basic physics as the friction resulting from 
the fluctuating gravitating forces acting on the satellite as 
a consequence of the non-continuos character of the particle 
system. The fact that fluctuating forces cause dissipation is 
quite a general scheme in Physics. It is at the basis of phys- 
ical phenomena such as electric resistance in conductors or 
viscous friction in liquids (Reif 1965). 

Dynamical friction has important consequences in the 
evolution of astronomical systems, mainly because it causes 
a decay of orbiting bodies, so that merger timescales and dis- 
sipation rates by dynamical friction are closely related. The 
merger scenario is at the basis of a great deal of processes 
in Astronomy. Not only is it the framework for the general 
galaxy formation picture in hierarchical cosmological mod- 
els, but also for more particular aspects of the evolution of 
a number of astronomical systems, such as galactic nuclei, 
cD galaxies in rich galaxy clusters, compact galaxy groups 
and so on. 

A dynamical friction formula was first obtained, in a 
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kinematical approach, by Chandrasekhar (1943). He has cal- 
culated the rate of momentum exchange between the test 
and field particles as the result of a sum of uncorrelated 
two-body encounters, obtaining the expression; 



M—^ . InAp < v)v, 

at v' 



(1) 



where M and v are the test particle mass and speed, re- 
spectively, p(< v) is the density of background particles 
with velocity less than v, G is the gravitation constant and 
A = , with Pmax and Pmin the maximum and minimum 
impact parameters contributing to the drag. Keplerian or- 
bits for both the test and field particles and an infinite and 
homogeneous background have been assumed in the deriva- 
tion of eq. (^) 

Chandrasekhar's formula is widely employed to quantify 
dynamical friction in a variety of situations, even if in most 
astronomical problems the background is neither infinite nor 
homogeneous. This formula is known to give the correct or- 
der of magnitude, but it suffers from several drawbacks, that 
arise from the very physical assumptions made in its deriva- 
tion. Furthermore, it cannot describe some situations, as for 
example the drag experienced by a satellite placed outside 
the edge of a finite gravitating system. In fact, the satel- 
lite would be decelerated on physical grounds, because it 
causes a perturbation to the system. However, according to 
eq. (|l|), the drag would vanish. For this reason, other works 
on dynamical friction followed Chandrasekhar's pioneering 
study, either from a numerical (Lin & Tremaine 1983; White 
1983; Bontekoe & van Albada 1987; Zaritsky & White 1988) 
or an analytical point of view. Analytical descriptions have 
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the advantage that they help understanding the underly- 
ing physics. Several methods have been developed: Fokker- 
Planck equation based (Rosenblunth, Mc Donald & Judd 
1957; Binney & Tremaine 1987) polarization cloud approach 
(Marochnik 1968; Kalnajs 1972; Binney & Tremaine 1987; 
Weinberg 1989; Bekenstein & Zamir 1991), resonant parti- 
cle interactions (Lynden-Bell & Kalnajs 1972; Tremaine & 
Weinberg 1984; Palmer & Papaloizou 1985; Weinberg 1986). 

Berkenstein and Maoz (1992, hereafter BM92) and 
Maoz (1993, hereafter M93) introduced a fluctuation- 
dissipation approach to dynamical friction. The fluctuation- 
dissipation theorem (Kubo 1959) relates the friction coeffi- 
cient to the time integral of the correlation function for the 
fluctuations causing the friction. They showed that dynam- 
ical friction fits into this general scheme, which provides a 
powerful technique to study it. This approach is in fact a 
return to Chandrasekhar's original attempt to give a sta- 
tistical description of dynamical friction. Other stochastic 
approaches to dynamical friction are those of Cohen (1975) 
and Kandrup (1980). M93 derived a formula for the drag 
experienced by an object which travels in an arbitrary mass 
density field, assumed to be stationary and formed by par- 
ticles much lighter than the object. 

A common feature of all the previous approaches is that 
they treat the satellite as being rigid and without struc- 
ture in the velocity space. This approximation can be good 
enough in a number of astronomical situations, but in oth- 
ers, these two ingredients could play a crucial role. As a first 
example of such a situation, let us consider the dynamical 
evolution of compact groups of galaxies (e.g. Mamon 1993) . 
In this case, the velocity dispersion of individual galaxies is 
comparable to the velocity dispersion of the common halo 
that hosts them. As a second example, we recall that in the 
problem of the interaction of two comparable mass galax- 
ies, the energy exchange due to dynamical friction cannot 
be calculated in the previous frameworks, because their pre- 
sumably comparable velocity dispersions need to be taken 
into account. 

In this paper we present a fluctuation-dissipation study 
of the orbital changes experienced by a non-rigid satellite, 
composed of gravitating particles with a finite velocity dis- 
persion, as it interacts with a general background. An exten- 
sion of BM92 and M93 techniques has allowed us to calculate 
the rate of energy exchange between them as a result of fiuc- 
tuations in the gravitational forces of both the background 
and the satellite. 

The paper is organized as follows: in Sect. 2, the phys- 
ical formulation of the method and the general expressions 
giving the energy exchange rate are presented. In Sect. 3, 
we calculate the instantaneous energy variations for general 
backgrounds at rest and with a Maxwellian velocity distribu- 
tion. Some particular limits are dealt with in Sect. 4. Finally, 
in Sect. 5, we summarize and discuss our results. Two Ap- 
pendices follow, where the results of the calculation of the 
correlation matrix and of an integral needed are given. 



2 PHYSICAL FORMULATION 

We will study the energy exchange between two self- 
gravitating equilibrium systems called the background (B) 
and the satellite (S). The background consists of Nb equal 



mass particles, each of them with a mass ms. Its total mass 
and typical size are Mb = itibNb and Rb, respectively. 
These particles exhibit a velocity spectrum with dispersion 
cr_B and zero mean at t = (which is equivalent to consid- 
ering the origin of the reference system placed at the halo 
center of mass at t = 0). The satellite in turn is composed of 
Ns equal mass particles of mass ms, total mass Ms, typical 
size Rs, and a velocity distribution with dispersion erg and 
mean equal to the center of mass velocity of the satellite, 

VCMS- 

As we are mainly interested in the effects that a non 
vanishing as would have on the evolution of the whole sys- 
tem, we can consider that as ^ ctb = a. The virial theorem 
tells us that a ~ GNmR~^ , so that as ~ aB implies that 

NsmsRs^ ~ NBrnBRB^^. (2) 

As both the satellite and the background are assumed 
to have a finite size, we can safely consider that 

Ni/-->(^^y (3) 

for whatever Q, /3, 7 = B or S, and n, m of order unity. If 
nia 3> rria (3 is the particle class contrary to a), as could 
be the case, for example, when the a and a systems are 
composed of stars and dark matter particles, respectively, 
then eqs. (^) and (^ imply the inequality Na <C Na- But, 
in any case, both Na and Na will be assumed to be very 
large. 

2.1 The fluctuating forces acting on a particle 

The fluctuating force, Fi^ (r^g , t), acting at time f on a satel- 
lite particle is placed at position rig and caused by its in- 
teractions with both background and satellite particles, can 
be written as: 

F,, (r,g , t) = Ffg (r,g , t) + Ff^ (r.^ , t) (4) 
where 




are the fluctuating forces acting on is caused by the back- 
ground and satellite, respectively. Each of these forces re- 
sults from subtracting to the total many-body or discrete 
force a smooth part derived from the mean field potentials, 
^^(r) and $s(r), due to the smooth densities, nB{r) and 
ns{r), respectively. These satisfy the relations: 

na{r) = Na I dufo{r,u), (7) 
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V^$q(t') = 4TTGmana{r) 
and 

drna(r) = N^, 



(8) 



(9) 



where a = B,S, /o'(r,ii) is the one-particle distribution 
function for background and satellite particles in the unper- 
turbed state, and the ^^r~^ factor takes into account that 
the is particle does not interact with itself at time t. An 
exchange of the B and S labels in eq. ^ gives the expres- 
sion for the fluctuating force on the background particle is 
caused by the satellite. Changing the S labels into B in eq. 
(P) we get the force on is caused by the background fluctua- 
tions. In a compact formulation, we can write the fluctuating 
force on a generic class a particle, ia, due to class /3 particles 
as; 



(10) 




where now Af^ — N, 



takes into account the possibility 
that a — p. 

The equation of motion of either a background or a 
satellite particle at time t can be written in a compact no- 
tation as: 



Jo 



(11) 



where Vi^.o is the particle velocity at time i = 0. 

Stochastic forces are weak as compared with the smooth 
global forces. In fact, the fluctuating force between a class 
a and a class l3 particle is of order (Kandrup 1980, Saslaw 
1985) 

\K,fiuc\=O{Gm„m0n^^'), (12) 

where n ~ N^Ra^ + (1 ^ Sap)Ni3R'^^ is the average particle 
number density for the whole system. The smooth density 
of class 5 particles causes a gravitating force on a class 7 
particle that is of order: 



O (Gm-fmsNsRs) . 
Eqs. ( |l^ and ( |l3| ) give, after some algebra: 



(13) 



= O 



\ms J I -R/3 



(l-5a 



\ms J 



3/2 



-1 2/3 



,(14) 



where eq. (g) has been taken into account. Now, if ms 
rriB-, then we wiU have 



let ■,fluc 



j-y , smooth 



= O 



|i) +(1-...) 



2/3 



-1/3 



.(15) 



this is indeed a very small quantity. If, on the contrary, back- 
ground and satellite particles have very different masses, eq. 
(jiM must be examined more carefully. The most unfavor- 
able case is when ms <^ m^,"f — S and a — [3 — 5. In this 
case, eq. ( p^ ) gives: 



1 , / Z u c 



= O 



4/3 



Rs_ 

Rp 



4/3 



-1/3 



and it would be sufficient to ask that 



iV 



1/3 



> 



(I 



\ms 



4/3 



(16) 



(17) 



to ensure that the fluctuating forces are much smaller than 
the smooth ones. 

Next we compare the timescale for fluctuation. Tap, and 
the timescale T^s over which the velocity of class 7 parti- 
cles changes appreciably as a result of the smooth forces 
produced by class 5 particles. The first is of order of the 
time required for the nearest neighbor to travel an inter- 
where a IS (J OL or (J a and 



particle distance. 



TaP = 



' Sai3)ni3. As it may happen that Ha and np are 
very different, we take the maximum Ta/s timescale, corre- 
sponding to = Nm,inRS'i„. The second timescale is set 

\V^l\Rs 



by = ^ 



SO that we can write 



T^S 



V Rs ) ' 



(18) 



where we have used the inequality (^). 

The rate of energy exchange between subsystems B and 
S will be obtained as an integral on time which involves 
the correlation matrix (see below and BM92 and M93). The 
correlation matrix is known to fall with time faster than 

(Cohen 1975; Kandrup 1980; M93), so that it can be 
taken to vanish for times much larger than the fluctuation 
timescale, and in particular for times of the order of the 
macroscopic timescale, T^s- As a consequence, in this work 
we will consider time intervals after t = 0, St, which are 
large as compared with the fluctuation timescale, t™^, but 
much shorter than the macroscopic timescale (see Reif 1965 
and M93 for a detailed discussion) 



'a/3 



(19) 



2.2 The energy exchange 

The energy of a particle ia is not conserved. Its total in- 
stantaneous variation at time t, due to the fluctuating forces 
caused by class P particles is given by: 



dt 



(20) 



Taking into account the expression for Vi^ (t) given by 
eq. (0), we get for times t verifying ( |l9| ) 



\ dt 



-Flit)- 



,0 + — > J dt'Fjit') 



-(21) 
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where the potential terms have been neglected because 
I WicO 1^1 V"I>7 1 t, as ensured by ([l9|). The total energy 
change of the ia particle during the time interval (0, 5t) is 
easily calculated by integrating eq. . Summing up on ia 
we get the total energy variation of class a particles due 
to the fluctuating forces caused by (3 particles in this time 
interval: 



pot 

+ —yfdt'Fj{t' 



(22) 



We will be interested in the average value of this quan- 



tity. 



2.3 Statistical averaging 

The presence and motion of the 15 satellite particle perturbs 
the background. As a result, while the statistical (ensemble) 
average of the background fluctuating forces acting on the 
satellite vanishes in the unperturbed state, {Ff^(t))o = 0, 
they do not vanish anymore in the real perturbed system, 
(Ff^{t)) 7^ 0. The same is true in general for the stochas- 
tic forces caused by class /3 particles and acting on class 
Q particles: (Ff^(t))o = 0, but {ff^(f)) / 0, because the 
phase density and the number of microstates available to 
the system has changed due to the perturbation. Following 
BM92 and M93, we assume that the probability of finding a 
dynamical variable, Q, with a given value, is propor- 

tional to the change in the number of microstates available 
to the whole system. This change is given by the factor: 



K = exp [SS] 



(23) 



where 5S is the entropy change of the system as a result of 
the distortion. The expected value of a dynamical function, 
Q, can now be written as: 



(Q) 



J drQ{r)fo{r)K 
J drfoir)K ' 



(24) 



where F are the variables defining the phase space of the 
Nb + Ns particles and /o(F) is the distribution function for 
the unperturbed state. 

As a result of the fluctuations, a generic particle initially 
with energy Si^ , will change to an energy e'i^ = + Ssi^ , 
a — B,S. Recalling the definition of entropy: 



S = So- / ,fln/dr. 



(25) 



with So a constant, this energy change results into a total 
entropy variation given by: 



(26) 



where /"(ei^) is the class a one particle distribution func- 
tion for the unperturbed state corresponding to an en- 
ergy Ei^ . Most two-body encounters are weak, and then 
Ssi^ ^ £,;„. Expanding the r.h.s. of eq.(E6|) and then the 



exponential in eq.(|23[), we obtain that the change in the 
probability function is given approximately by: 

1 dr 



K 



EE 



(27) 



For most astronomical applications it is justified to use 
an isothermal Maxwellian one-particle distribution function, 

where u is the velocity vector, aa is the velocity dispersion of 
class a particles and f3a = ^ ''"^2 is an inverse temperature. 
Introducing this expression for the distribution function in 
eq. ([2 7]), the factor of change becomes: 



(29) 



The next step is to find out an expression for the energy 
variation. The total energy of the system must be conserved, 
so that we have: 



EE 



= 0. 



(30) 



Moreover, for one given particle, ia, its energy variation 
can be expressed as: 



(31) 



where ai^ is the fraction of the total class a autointeraction 
energy, AE^.tot, absorbed by particle ia, and hi^ is the frac- 
tion of the total interaction energy between class ct and a 
particles, caused by the fiuctuating forces of the a subsys- 
tem, absorbed by ia particle. Summing on ia in eq. (^l|), we 
get the total energy change of subsystem a: 



E 



5£ic = ^Ea,tot ~ ^E^^tat, 



(32) 



compatible with energy conservation (eq. (|3C|)). To write 
down eq. ( ^2[ ) from (pl|), it has been taken into account 
that = bi^ — 1. Energy signs are such that 

AE^' > if particle ia gains energy due to the fiuc- 
tuating forces caused by class 7 particles and conversely. 
With this convention, if energy fiux is from subsystem a 
to subsystem a, then AEa^tot < (subsystem a loses en- 
ergy), and AE^ t^t > (subsystem a gains energy), so that 
y^^- Ssi^ < and fei- > 0. Inserting eq. ( p2[ ) in eq. 
( [29I ) and recalling eq. (^2|), we obtain: 



- i;i:/'-«n(.')]f 



(33) 



Eq. ( |24| ) allows us now to write the ensemble average 
of the instantaneous energy variation at time t of a generic 
class fj, particle due to the stochastic forces caused by class 
V particles: 
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V 

dt 



dr/o(r) 




(34) 



To second order in the fluctuating forces, taking into 
account tliat the average of the stochastic forces in the un- 
perturbed state vanishes, and summing on the subindex 
to obtain the global effect, we get: 



dt 



E- 



dsTr 



, J -St 



where 

Cr!,As)^{FUQ)®FUs)), 



(35) 



(36) 



is the correlation matrix (see Appendix A; the symbol (8) 
stands for the tensorial product of vectors f °^ and -F"^^, 
the average is on the unperturbed states of particles a and 
/3 and the invariance of the correlation matrix under time 
translation has been taken into account). Because in the 
unperturbed state two different particles are not correlated, 
the correlation matrix defined in the previous equation van- 
ishes if a 7^ /3. This has been taken into account to deduce 
eq. (^5|) from eq. (p^. Note that because for s > 5t the 
correlations vanish, the lower limit of the integrals can be 
extended up to — oo. Eq. ( ^sj) is the expression of the global 
instantaneous energy variation of subsystem /i due to the 
fluctuating forces of subsystem v. The integrand in the first 
term of the r.h.s. of eq. ( ^sj ) is invariant under time rever- 
sal. Extending the integral to positive s, this term is the 
corresponding power spectrum at zero frequency (Wiener- 
Khintchine theorem, see Reif 1965 and BM92) and, conse- 
quently, it is positive and represents a heating term of class 
H particles due to the fiuctuating forces caused by u parti- 
cles. It is of order 0{-^) ^ 1 relative to the second term. 
Terms of this kind will be neglected in this work. The sec- 
ond term can be either positive or negative, depending on 
the sign of (/?£; — Pv). In the next section it will be explicitly 
calculated for = B, S. 



3 THE INSTANTANEOUS ENERGY 
VARIATION 

3.1 Energy exchange between the satellite and 
the background 

Energy fiows between the satellite and the background as 
a result of either the global energy variation of the satellite 
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particles due to the fluctuating forces of the background, or, 
conversely, the global energy variation of the background 
particles due to the fluctuating forces of the satellite. These 
fluxes are given by eq. (^Hj) when ^ ^ v. 

Let us first calculate the effect due to the stochastic 
forces of the background. This term could cause an energy 
flux responsible for the satellite deacceleration. Eq. ( |35[ ) with 
H = SyU = B and neglecting the heating term, gives: 



\ dt 



= {Ps-Pb) 



EE 

's is 



dsVisfiCf^%{s)vjsfi, 



(37) 



where the correlation matrix is given by eq. ( AlO ) with a = 
B and 7 = 5 = S. In the case of a point-like satellite with 
ms ^ TUB, Ps Pb and then the satellite loses energy to 
the background. In our case, however, it is not excluded in 
principle that Ps > Pb and then the energy would flow from 
the background to the satellite. 

The integral over ds in eq. (^^ can be calculated taking 
into account eq. (AlC) and the equality: 



ds; 



A + Vs 
A + Vs\ 



A-V + AV 



A 

1 



V 



(38) 



where A± is the projection of vector A on the plane normal 
to vector V . In our case either V = Vj„ — v^g or V = Vj„ , 



(see eq. ( |A10[ )), and A 
obtain: 



or A 



We 



\ dt 



r~i2 2 2 , r, a \ ^B 

= G rriBmsiPs - Pb)j^ 



EE 



E 



is is L h 
■"is ■ (r 



3s -rheWr- 



2 



'■js - 11 «js - ?\ {vjs - rhs)± \ 
rjs -rhg P (wjs -«/iB)-'"is 



^JS '"'hi 



V 



3S 



1^1 (rjs -rhB)± P 



+ 



■Thf 



'"hB 11 {r^s -rhB)± 1^ 



/ drdr' nB(r)nBir')- — ^^^'^ 
Nb. ri„ - r 



(39) 



where the sum on a generic subindex, ia, means Na times 
the average on velocity and positions of class a particles, and 
is carried out by means of the distribution function given in 
eq. (g). 

To proceed further, we recall that the velocity distri- 
bution of background particles is assumed to be Maxwellian 
with zero mean and dispersion as, and the velocity distribu- 
tion of the satellite particles has mean equal to the satellite 
center of mass velocity, vcms, and dispersion os- The inte- 
gration over dvhB gives: 



dt 



G mgrnsiPs ~ Pb) 
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EE 



drhgUBirhB) 

is is 

exp(7fs - y|g)erf(7j.;) - 1 
I ^is ^ ^hf? 



I ^1 I drdr'nB(r)nB{r')- ' 



Nb 



ris 



r< 



(40) 



where erf stands for the error function, y. 



x/2<y 



lis 



^''"^ and taking = 1. 



Vjs ■ C with C = |r,-^-r,^| — ~^ 

The integration over dvj^ involves the integral: 



In Appendix B we show that 

2 

/ = -f- exp(QcA/s - xcMs)eri{acMs), 



(41) 



(42) 



where = (t|+o-|, xcms = and acMS = xcms'C- 

The integration over dvi^ is trivial recalling that Vi^ — 
vcMS + and that for a Maxwellian distribution the 
contribution of Ui^ to the first moments vanishes. Taking 
maria = Pq, we finally get: 



dt 



-fUl^^^lfdr dr 



drhsPB{rhs) 



jsPs{ris)psirjs) 
(ris - Vhs) ■ vcMS 

I f^i Q hfi I 



X <^ -^exp{acMS - s;cMs)erf(acMs) - 1 



1 / , , / / N / /N [{ris - r) ■ vcMs] 
+ -r^ / drdr pB{r)nB[r )- 

Nb 



- r 



(43) 



A specification of the density distribution of both the 
satellite and the background is needed in order to carry out 
the integrals over the space variables. 

Next we calculate the variation of the background en- 
ergy caused by the stochastic forces of the satellite. It can 
be obtained from eq.(p5|) with fi — B and ly — S, neglecting 
the heating term: 



dt 

= (/3b ~Ps)Y,Y. dsv,BS>c!^% {s)vjB,o. (44) 

iB is •'-'^ 

For a background at rest with a Maxwellian velocity 
distribution function, this energy rate vanishes when one 
performs the integration over dvi^ . We conclude that the 
effect of fluctuating forces of the satellite acting on the back- 
ground only heat it, and have no effect on a variation of the 
orbital energy of either the satellite or the background. 

The total instantaneous energy flow between the satel- 
lite and the background is given by the difference between 
the l.h.s. of the eqs. (||) and (Q) (see eq. (||)). It can be 
written in terms of the rate of change of the satellite orbital 
and internal energies as; 



dt 



dt 



+ 



f dEg 



\ dt 



dE 



S.tot 



dt 



(45) 



where -Ef^b ^i^d Ef^ stand for the orbital and internal energy 
of the satellite, respectively, and the second equality results 
from the zero value of the flow given by eq. (^). The rate 
of change of the satellite internal energy can be calculated 

with the help of eq. (^ij) with Q — — or: 



(46) 



where Ui^ = Vi^ — vqms is the velocity of is particles with 
respect to its center of mass. This gives an expression similar 
to eq. (k3), except that now vcms takes a zero value, so that 



I dEf„ 
\ dt 



= 0, 



(47) 



and then eq. ( }43D gives, at second order in the fluctuations, 
the rate of change of the satellite orbital energy. 

3.2 The self-interaction energies of the 
background and the satellite 

When fj, = 1/ — B OT S, eq. ( |35[ ) describes the instantaneous 
energy variation rate of class /i particles due to the stochastic 
forces caused by ^ particles themselves. 

According to the eq. ( |3^ these energies do not play any 
role in the energy change of the subsystems S or B. They 
only represent the energy change of an individual particle 
(see eq. (|l|)). 

The self-interaction energy of the background is easily 
obtained from eq. (pq) with fj, = u = B: 



dt 

= (/3s - /3s) E E / dsv,^,oCff,^{s)v,„o. (48) 



Again, the integration over dvi^ makes it vanish for 
a background at rest with an isotropic velocity distribu- 
tion function: because no translation energy is available, the 
autointeraction results only in a slow heating of the back- 
ground. 

Regarding the autointeraction energy of the satellite, 
eq. (pq) gives: 



dEn 



dt 



{13b - l3s)Y,Yl / dsv,,,oC! 

'i —oo 



ss 
siE 



"iB.o- 



(49) 



The correlation matrix is given by eq. ( AlO ) with 
a = Sj'y — S and 5 — B. Substituting the expression for 
Cffjg (s) in eq. (^^ , and performing the integral over ds 
with the help of eq. (pSl), we obtain: 
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dt J^J Ns 



EE 



E 



(T-jg - Ths) ■ Vis 
I '''is ^ ''hs I'' 



is JB L hs 

\riB - r,,^ 11 - Vhs 1^1 (r-jB - r-h.Jx | 

I T-jB ^T-fas P Kb -"fcs)-«jB 

'■jB -■'■hs II '"iB -fhs PI (t-jb -rhs)± P 

I'jB -«fes 11 (rjB -~rhs)± P 
+ Tr / '^^'^^ ns(r ns(r )- — 



.(50) 



The integrations over velocities can be carried out fol- 
lowing the same steps as in the previous section and we 
finally have: 



dEs,tot\ \_fmso-g \ G 



dt J^,/ \rnBO-g 



drisdrjgPsiris)pB{rj 



-4 fdr,,ps{r,s)^^"''~"''^ 



The 



/ drdr ps{r)ns{r ) ^ 



rrs - r PI »-jB - r' 



(51) 



where acMS = xcms ■ ,r -r \ ■ 



4 PARTICULAR LIMITS 
4.1 Point-mass satellites 

In eq. (^) it is implicitly assumed that the distance between 
a generic satellite particle and a generic background parti- 
cle cannot be smaller than a scale, dmin, that is, that there 
exists a minimum effective impact parameter. As the main 
contribution to the integrals appearing in this expression 
comes from small j rg — rs 1 values {rs and rs are generic 
satellite and background particle positions), an accurate de- 
termination of the dmin scale is a crucial point when study- 
ing dynamical friction. This scale, however, cannot be deter- 
mined by the present approach to this problem. There exist 
in literature several estimates of d,nin for specific situations 
(White 1976; Bontekoe & van Albada 1987). White (1976) 
in fact shows that for spherical symmetric satellites, it is a 
good approximation to consider them as point-like systems 
with a cut-off in their distances to background particles, 
and that for the case of satellites with a King mass distribu- 
tion (King 1966), this cut-off is given by dmin — where 
Rt is the satellite tidal radius. A point mass satellite has a 
density profile given by ps{rs) = MsS{rs — rcMs), where 
5{rs — rcMs) is the delta function in three dimensions and 
rcMS is the position of the satellite center-of-mass. With 
this value of psifs), eq. (fl3|) gives in the point mass limit: 



dt 



drhBpB{rhB)'S>{rhB) 



{rcMS - rhB ) ■ vcMS 



rcMS - 



^2 

-^exp{acMs - a;cMs)erf(acMs) - 1 

-|-— |— / drdr' pB{r)nB{r')Q{r)Q{r') 
Nb J 

[{rcMS - r) ■ vcMs] 



rcMS 



rcMS 



(52) 



where now acMS = xcms ■ il*^"'^^ ^''^i s^nd 6(r) 

' CMS — ' Hb I 

6 (I rcMS — r \ —dmin) is the step function. 



4.2 Homogeneous backgrounds. The 
Chandrasekhar limit 

When the background density is uniform, we can set 
Ps(»'b) — pB,o and rcMS = in eq. (^) and then it be- 
comes: 



orh 



X [eri{xcMs) 



mBcy% 
msiy% 



- 1 



47rG^M|ps,olnA 



VCMS 



^XCMS e:^p{-XcMs)\- 



(53) 



This recovers the Chandrasekhar formula for the motion 
of a massive test particle in an homogeneous background, ex- 
cept for the factor containing the temperature ratio (which 
vanishes in the Chandrasekhar formula because ms ^ tub 
in this limit), and, now, or = (o-| + cr^)^^^ at the place of 
the background velocity dispersion. 

In the Figure |l| we plot the ratio, R, of satellite en- 
ergy loss in the Chandrasekhar limit with as = and with 
as 7^ 0, assuming that Tb/Ts = 0. As can be seen in this 
Figure, the effect of a non zero as increases with increasing 
ITS and is more important at low vcms- The Chandrasekhar 
formula always overestimates the dynamical friction force, 
and in some situations the effect of neglecting the satellite 
velocity dispersion could cause an error as high as an order 
of magnitude. 



5 SUMMARY AND DISCUSSION 

We have derived an expression giving the orbital energy 
exchange due to dynamical friction experienced by an ex- 
tended body, composed of Ns bound particles endowed with 
a velocity spectrum, as it moves interacting with a non ho- 
mogeneous discrete background. It has been assumed that 
both, the satellite and the background, have Maxwellian ve- 
locity distributions and that the background is static. 

Self-interactions of both satellite and background par- 
ticles have been taken into account. This results in no effect 
on their energy exchange. 

Heating terms appear in quite a natural way in our 
approach both due to interactions among particles of the 
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0.5 1 1.5 Z 2.5 



Figure 1. Ratio, R, of satellite energy loss in the Chandrasekhar 
limit with 0-3 = and with erg 7^ 0, assuming that Tg/Tg = 0. 
Lines correspond to different us I as ratios from 0.2 to 2.0 (ratio 
increases as the thickness of the line increses) 



same kind or of different kind. They are a factor of XjN 
smaller than the orbital effects. 

In the point-like satellite limit (or small as compared 
with the background size) and constant density background, 
we obtain an expression that recovers, in the limit as/uB 
and Ts/Tb 00, Chandrasekhar's dynamical friction for- 
mula (eq. (|l|)). Its comparison with the energy loss given by 
eq. (|l|) allows for a quantification of the effects of having 
a non zero as- It has been found out that the energy loss 
is always smaller in this case, and that the difference can 
be up to about an order of magnitude for slow satellites as 
compared with crs. 

In deriving eq. (js^) we have considered time intervals, 
St, that are short as compared with the time scale for the 
variation of the particle velocities and positions. This al- 
lows us to neglect the effects of the smooth gravitational 
potential gradients. However, in order to carry out a pre- 
cise calculation of dynamical friction, the whole history of 
the system from an early enough time and the interactions 
along the entire satellite trajectory should have been taken 
into account. This would have made the problem extremely 
difficult to solve. Instead, taking only a finite St, means that 
interactions with distant particles have not been accurately 
considered. Nevertheless, we recall that the contribution of 
particles to dynamical friction quickly decreases with dis- 
tance, so that this neglect should not result in major conse- 
quences. 

The bound of St has also another consequence: this ap- 
proach is unable to describe slowly accumulating effects on 
dynamical friction (Kalnajs 1972) or the effect of reversible 
dynamical feedback (e.g. Tremaine & Weinberg 1984), be- 
cause they arise as a consequence of the periodic motion of 
the satellite after many revolutions. 

Despite these shortcomings, the extension of the 
fluctuation-dissipation approach to dynamical friction pre- 



sented in this paper, has resulted in the derivation of a for- 
mula that takes into account the space and velocity struc- 
ture of the satellite. This represents a common situation in 
many astrophysical processes and, as we have shown, might 
have important quantitative consequences in the setting-up 
of timescales for these processes. 
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APPENDIX A: THE CORRELATION MATRIX 

The correlation matrix is defined in eq. ( ^^ with the fluc- 
tuating forces given by eq. (|lo|). Each matrix is an object 
with four index, and each index takes on two different val- 
ues {B and S). This makes sixteen different possibilities, 
corresponding to the tensorial product of four different pos- 
sibilities for the Ff_^ (t) forces. Writing the Fourier transform 
of I r — r' and na{r) with respect to r, the convolution 
theorem and eq. (0) imply that: 
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J2 expHfc-r-,Jt)]-i^. 

where 

n'j^ = J drna{r)exp[—ik ■ r]. 

With this definition 
(exp[-ifc ■ riJ)o = 



NL 



(Al) 



(A2) 



(A3) 



and then the average of the stochastic forces, eq. (Al), in 
the unperturbed state vanishes, as required. 
Inserting eq. ( [Al[ ) in eq. (p^), we get: 



G mami3m^ms 



dk dl 



47r4 

k ® I exp[i(fc • Ti^ {Q)+l-rj^ (s))] 



/ , exp[-ife • r9„(0)] - ^rife 



^ exp[-ii • rhi, (s)l - -^nf 



(A4) 



and eq. (A3) now gives; 



47r4 



/ ^^fe®'exp[i(fc-r,^(0)+rr,,(s))] 



(A5) 



where 

Eki{ga,hf3,a) = {exp{-j[fe ■ (0) + / ■ rhj3(s)]})o. (A6) 

It is assumed that particles are uncorrelated in the un- 
perturbed stat e. Th en, if a 7^ /3, necessarily ga 7^ hp and the 
average in eq. (A6) factorizes, as corresponding to uncorre- 



lated variables, giving Cf^-^ (s) = 0. The average also fac- 
torizes for different particles belonging to the same particle 
class, i.e., when a = /? but ga ha- There are N^iN'^ — 1) 



such terms, and each of them has the value 



When 



get = ha, that is, when we consider the same particle, the 
average does not factorize anymore. Taking this considera- 
tions into account, the correlation matrix becomes: 



^2 2 



dk dl 



k ® I exp[i(fc ■ (0) + r (s))] 



N'aEu{ga = ha,s)-N, 
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(A7) 



The position of a generic particle at time s satisfying 
( p^ can be written as: 



(A8) 



An integration of eq. ( |ll[ ) would give four more terms 
corresponding to the gravitational acceleration and the fluc- 
tuating forces. The acceleration terms are negligible when 
compared with Vj^^s (see eq. (p^) and the fluctuating 
forces would give rise to third order terms. 

Once this expression for the particle positions at time s 
is substituted in eq. {^7), the integrals over dk and dl can 



be easily calculated taking the gradient of the Fourier rep- 
resentation for the Green's function of the Laplace equation 
in three dimensions: 



dk .2ttA 
— fcexp[jfc- A] = 1 — . 



(A9) 



Finally, combining eqs. (|A2|), (|A^), (|A7|), (|A^) and (|A9|) 
the correlation matrix reads: 



Caa I \ ^2 2 
J jAs) — G rriamryms 



K 

Na 



(r^ -rhj [rj, - r^^ + s{vj, - vhJ] 



J drdr'ua {r)na{r') 



r + svj 



r' + SDo 



(AlO) 



The subindex have been dropped from the r and v 
vectors. 



APPENDIX B: 

In this appendix we calculate the integral 



^ ^ / '^■"/m(«) exp [7^ - y^] erf(7), 

where v is the velocity of class S particles, that is 

V = VCMS + U, 



(Bl) 



(B2) 



whose distribution function, f^j, is a Maxwellian isotropic 
in u 

\2 " 



(27ro-|)3/2 



exp 



{v - VCM s) 

2ol 



(B3) 



In the last equation y = v/\/2aB, 7 = 2/ ■ C) and ^ is a 
unit constant vector. Expression (E3) implies: 



7 = 7ci\/s + lu 
where 

VCMS , 



(B4) 



7CMS 



and 



V2i 



UCMS A 

1-^ = — c 



^/2< 



The definition of the error function and equation (B4) 
imply that: 
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2 



erf(7) = erf(7cA/s) 



+ 4= / dtexp [-[t + 'YCMsf] ■ (B5) 
V Jo 

We now write / (eq. (Bl)) as an integral with respect 
to the relative velocities u. Taking ^ — (0, 0, 1) we obtain: 



I = 



exp (acMS-xcMs) 



(27r4)3/2 



dUx exp {Aux + Bx)^]^ 



erf(7cMs) + / dt exp [- {t + 7cms)^ 



(B6) 



where xcms = vcms olcms = Xcms ■ Cj ''"t = 

a% + cr|, A = ar/V^o-scrB, Bi = ^xcmsa, with i = x,y, 

and 7u = Uz/y/2TvaB- The integrals over du^; and duy can 
be easily performed and they give; 
^ n 2 



■erf(oo) 



TT 
l2 



The integral can now be written as: 
^2 |- 2 - 

TT 



erf(7CM5) + -/(a,fo) exp TacA/s - a^cAfs] (B7) 



where 



7(a, b) = — ds exp 



dtexp [-(t + (B8) 



with a = OS /oh, b = jcms, comes from the integration on 

The I{a, b) integral can be evaluated as follows: first, we 
derive it with respect to b, the we carry out the integration 
with respect to t, obtaining: 



dl{a,b) 
db 



= \/7r 



1 / b" 

(l + a2)i/2 1 i + „2 



— exp(— 6 ) 



(B9) 



and, finally, an integration on b leads to: 



/(a,fe) =/(a,0) + - 



erf 



(l + a2)V2 



erf(6) 



(BIO) 



with /(a, 0) = because erf(s) = — erf("s). Substituting the 
expression for I{a,b) in eq. (B7), we get: 



2 

/ = -f- exp (acAfS - a;cAfs) erf(Q:cAfs) 



(Bll) 



where now xcms = vcms/V^ot and acMS = xcms ■ C- 



